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, Some known solutions for static charged fluid spheres of Tolman-VI type in general relativity are reexamined. The 
gravitational mass that appears in the Pant and Sah [1] and Tikekar [2] solutions is shown to be of electromagnetic 
origin in the sense that the gravitational mass, along with all other physical quantities, depends on the electromag- 
netic field alone. The energy condition, singularity and stability problems of the models are discussed thoroughly. 
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1. Introduction 

.Studies of static solutions in connection with stellar in- 
teriors have been of continuing interest for many re- 
iSearchers in the framework of general relativity. These 
^vorks can be divided into two categories: (1) static 
fluid spheres without charge and (2) charged static fluid 
'spheres. Some remarkable works on the neutral cases of 
^Einstein's field equations, related to the Schwarzschild 
'solution, were done by Tolman [3], Wyman [4], Leibovitz 
[5], Whitman [6] and Bayin [7]. On the other hand, so- 
lutions of the coupled Einstein-Maxwell field equations 
have been obtained by several authors. Some of the 
interesting works in this direction with the Reissner- 
Nordstom space-time can be found out in the papers by 
Weyl [8], Majumdar [9], Papapetrou [10], Bonnor [11], 
Kyle and Martin [12] and Cooperstock and de la Cruz 
[13]. 

The reason behind the work on the coupled charge- 
matter distributions is mostly to avert the singularities 
that occur in the Schwarzschild-like solutions. These 
singularities are of two types: (1) the coordinate singu- 
larity at r = 2m which provides the event horizon of a 
black hole of mass m, and (2) the physical singularity 
at r = , the centre of the spherical system. Now, in 
connection with the singularity, it is observed that the 
gravitational collapse of a spherically symmetric distri- 
bution of matter to a point singularity may be avoided 
in the presence of charge. This happens since the re- 
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pulsive Coulomb force due to the charge (in addition to 
the thermal pressure gradient in the fluid) opposes the 
gravitational attraction. Therefore, the issue of study- 
ing the Einstein-Maxwell space-time has always received 
considerable attention from the researchers [1,14-16]. 

However, there was another issue in studying charged 
static fluid spheres, which is associated with electro- 
magnetic mass models. According to Lorentz [17], in 
an extended electron "there is no other, no 'true' or 
'material' mass" , and this provides only an "electro- 
magnetic masses of the electron". Wheeler [18] also 
believed that the electron has a 'mass without mass', 
while Feynman [19] termed this type of models the "elec- 
tromagnetic mass models" , where mass originates from 
the electromagnetic field alone and hence gives a phe- 
nomenological relationship between the gravitational 
and electromagnetic fields [13,20-24]. 

Therefore, motivated by the above ideas, we have 
recently [25] investigated the class of astrophysical so- 
lutions previously obtained by us [26] and also by Pant 
and Sah [1] for a static, spherically symmetric Einstein- 
Maxwell space-time. While discussing the Pant and Sah 
solutions [1], which represent the charged analogue of 
the Tolman VI solution [3], we [25] have studied the 
case n = only related to the electromagnetic mass 
models [17,19], where n is a free parameter appearing in 
the set of solutions. In the present investigations, how- 
ever, it is shown that the electromagnetic mass models 
can also be constructed not only for the specific choice 
n = but for the whole range of < n < 2, even 
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including fractional values. 

It is known that the Tikekar solutions [2] represent 
a generalization of the Pant and Sah solutions [1], so we 
have also studied them as a continuation of the verifica- 
tion scheme. It is possible to show, from the expression 
for the charge and the equation of state, that the grav- 
itational mass depends on the electric charge alone. 

In the concluding part, it is discussed that there are 
some other solutions available in the literature where, 
especially in the case (A,n) in section X of Ivanov's 
paper [16], our requirement for a solution to describe 
an electromagnetic mass model, is also satisfied. For 
physical validity related to all cases studied in this arti- 
cle, the singularity problem, stability situation and the 
energy conditions are discussed in connection with de 
Sitter space-time. 

2. The Einstein-Maxwell field 
equations 

Consider the spherically symmetric line element 

ds^ = Qijdx'dx^ = e-'^^'Ut'^ - e^^^'Ur'^ - r'^dil'^, (1) 

with dfl'^ = dO^ + sin^ 9 dcfP' in the standard coordinates 
2,(04,2,3) ^ 1^^^ 0^ ^ where v and A are two functions 
of the radial coordinate r . 

The Einstein-Maxwell field equations corresponding 
to spherically symmetric static charged source are then 
given by 

e-^{>! /r^l/r^) + l/r^ = %-Kp + E^, (2) 
e-^{v' /r + l/r^)-l/r^ = d,TTp- E"^, (3) 
e-^K/2 + i/'V4 - v'X'/A + {v' - \')/2r] 



87rp + 



A/2 



(4) 
(5) 



where p, p, E and a are the matter energy density, the 
fluid pressure, the electric field intensity and the electric 
charge density, respectively. Here, the prime denotes a 
derivative with respect to the the radial coordinate r . 

Eq. (5) can be equivalently, in terms of the electric 
charge q, expressed as 



q{r) = r^E{r 



) = / 47rr^cr( 
Jo 



'.^/^dr. 



(6) 



Therefore, Eq. (2), by the use of the above equation (6), 
reduces to 



e"^ = 1 - 2M{r)/r, 



(7) 



where M(r) is the active gravitational mass, which can 
be expressed in terms of the effective gravitational mass, 
m(r) , in the form 

m(r) = M(r) + n{r) = An [ [p + q"^ /Snr^jr'^dr, (8) 

Jo 



/i(r) being the mass equivalent of the electric field. It is 
to be mentioned here that customarily people consider 
/i(r) as the m,ass equivalent of electrom,agnetic field be- 
cause of the fact that it is associated with the idea of 
Lorentz's conjecture of 'Electromagnetic mass' for the 
extended electron [17,19]. However, the models consid- 
ered here being static, the fields arc not electromagnetic 
but only electric. 

It is interesting to note that Eq. (8) provides the in- 
crease of the total gravitational mass due to inclusion of 
the charge [13,20,27], and in the absence of an electric 
charge it reduces to the usual active gravitational mass. 
In this regard, it is also to be noted that Eq. (8) is in 
the form of Bekenstein's [28] Eq. (31), where the part 
of mass related to the fluid is irreducible and charge- 
independent. Therefore, in the next section, one of our 
aims is to find out the criteria for constructing electro- 
magnetic mass models by showing that the matter part 
depends on the charge, and thus the entire gravitational 
mass turns out to be completely electromagnetic by ori- 
gin. 

3. Gravitational masses of purely 
electromagnetic origin 

3.1. The Pant-Sah models 

For a static, spherically symmetric distribution of charged 
fiuid. Pant and Sah [1] obtained a class of solutions 
which are as follows: 

= &r2", (9) 
e"^ = c, (10) 
P = 7^li-c(n-ij'J, (11) 

(12) 

(13) 
(14) 

(15) 
-1 
(16) 

With this expression of c, the above class of solutions, 
for the cosmological constant A = and also for the 
integration constant i? = 0, represents a charged ana- 
logue of the Tolman VI [3] solution. It is also to be 
noted that the matter density, the fluid pressure, elec- 
tric charge density and the radial electric fleld involved 
in Eqs. (11)-(14) follow the inverse square laws. Conse- 
quently, these physical parameters increase as the dis- 
tance decreases and become infinite at the centre of the 



P = 



167rr^ 
1 

167rr^ 



477 j-2 



[l-c{n-l)\ 
[c(n + 1)2-1], 



-{l-c(l-|-2n 



n^)} 



1/2 



E^ 
where 
h = 



2r2 



[l-c{l + 2n-n% 



-2n 



2m 



+1 

~ 9 



2m 



1 - 



2q' 



{l + 2n-n^y 
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spherical system. Therefore, due to singularities at the 
centre, such models are not physical as possible descrip- 
tions of stellar structure. They may be cither regarded 
as the Reissner-Nordstom black hole solutions, or one 
may think of a technique by which the singularity can 
be averted. We shall consider this point later. 

Now, the total gravitational mass m{r = a) can be 
calculated from Eq. (8) using Eqs. (11) and (14): 

n{2 - n)a'^ + 2q^ 



m ■■ 



2(1 + 2n - n2 



(17) 



A close observation of Eq. (17) shows that the status 
of the gravitational mass depends on the parameters 
n, a and q. Therefore, once we fix the values of a 
and q for a given charged spherical system, n can be 
regarded as an adjustable parameter 'having real, not 
necessarily integral, values' [3]. However, for physical 
viability, the values to be assigned to this parameter are 
< n < 2 . One can obviously assign innumerable values 
to n within this range. Related to the case n = , it 
has been already shown by us [25] that for vanishing 
electric charge all the physical quantities, including the 
gravitational mass, vanish, and the space-time becomes 
flat. We therefore, for our purpose, would like to study 
the following four subcases A-D only. 



A. 



0.5 



This case is favourable for some historical reasons. This 
choice was originally made by Tolman [3] himself in his 
uncharged version. In this case, the gravitational mass 
becomes 

3a2 + 8q^ 



m 



14a 



(18) 



The functional condition containing the energy density 
p and the pressure p, by virtue of Eqs. (11) and (12), 
can be written as 

n(a2 - 2g2) 



p + p 



1 



(1 + 2n - n2)a2 



which for the present case reduces to 



p + p = 



B. n = 1 



1 - 



2g2 



(19) 



(20) 



For this choice, the gravitational mass becomes 



m = 



2g2 



4a 



(21) 



whereas the functional condition containing the energy 
density p and the pressure p is 



p+p = 



87rr2 



1 - 



(22) 



It has been pointed out by Herrera and Ponce de Leon 
[27] that this particular solution, n = 1 , admits a one- 
parameter group of conformal motions and is homoth- 
etic. 



C. n= 1.5 

In this case, the gravitational mass is 
3a2 + 8g2 



m : 



14a 



(23) 



and the related functional condition here becomes 



P + P 



147rr2 



1 



(24) 



The gravitational mass in this case is the same as for 
n = 0.5. 

D. n = 2 

The gravitational mass in this case becomes 



r 

a 



(25) 



This expression is precisely the same as that in the n = 
case of our previous work [25] , it vanishes for vanishing 
electric charge and thus provides an 'electromagnetic 
mass' model [17,19]. However, the functional condition 
for the present situation differs from that of the n = 
case and is given by 



P + P- 



1 



27rr2 



1 



(26) 



Let us now, following Ray and Das [25, 26], and 
Ivanov [16] , choose a relation between the electric charge 
and the radial coordinate of the fluid distribution as fol- 
lows: 



q[r) = Kr'^ 



(27) 



where K and s are two free parameters. For the specific 
choice K = 1/^/2 and s = 1, the ansatz expressed in 
Eq. (27) reduces to q{a)/a = l/-\/2, where a is the ra- 
dius of the charged sphere. It is interesting to note that 
for this charge-radius ratio all the above perfect fluid 
functional conditions reduce to the form p+p = . This 
is known as the 'pure charge condition' [22] and also the 
imperfect-fluid equation of state in the literature since 
the matter distribution under consideration is in tension 
and hence the matter is named as a 'false vacuum' or 'de- 
generate vacuum' or 'p-vacuum' [29-32]. We would like 
to comment here that this imperfect-fluid equation of 
state precisely describes the cosmological constant, viz., 
the vacuum energy, and by the phrase 'pure charge con- 
dition' Gautreau [22] indicates the necessary condition 
for the 'Lorentz-type pure-charge extended electron' to 
be of electromagnetic origin. 

The above charge-radius ratio in turn makes all the 
total effective gravitational mass vanish for vanishing 
charge. Therefore, the gravitational masses here are of 
purely electromagnetic origin. It is to be noted here 
that all the above cases n = 0.5,1,1.5, including our 
previous case n = [25], provide 'electromagnetic mass' 
models with the imperfect-fluid equation of state. The 
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only exception here is the last case n = 2 where we 
have obtained an electromagnetic mass model even un- 
der the perfect fluid condition. Although for the charge- 
radius ratio q{a)/a = l/V^, the perfect fluid functional 
condition in Eq. (26) reduces to the imperfect one, but 
this ansatz is not necessary for making the gravitational 
mass in Eq. (25) suitably vanish as in the other cases. 

Let us look at Eq. (26) differently. The components 
of this equation, viz., the matter energy density and the 
fluid pressure are given by 



1 



P = 



P 



1 



87rr2 



(28) 
(29) 



It can be checked that for a positive pressure the con- 
straint on the charge-radius ratio to be imposed here 
is q{a)/a < ±2/3. The pressure vanishes for the value 
±2/3, whereas it becomes positive for values less than 
±2/3. Thus we have come across a case for electro- 
magnetic mass models where the equation of state is a 
perfect-fluid one with a positive pressure. This case is 
obviously in contradiction to Ivanov's observation [16] 
that "...electromagnetic mass models all seem to have 
negative pressure". This particular aspect has been 
pointed out in our previous work [25] where we have 
mentioned that there are some examples of electromag- 
netic mass models where positive pressures are also 
available, to be shown elsewhere. However, a few more 
examples in this direction need a further study. 

3.2. The Tikekar models 

We are now interested in discussing the Tikekar [2] so- 
lutions which represent a generalization of those of Pant 
and Sah [1]. As a particular solution, which describes 
a physically plausible distribution of a charged perfect 
gas, Tikekar [2] has obtained the matter density and 
fluid pressure in the following forms: 



(30) 
(31) 

(32) 



47r[(6 + d+l)2-4d]r2' 

d 



P 



47r[(6 + d + 1)2 - 4(i]r2 ' 
and hence the functional condition is given by 

b + d 



P + P ■■ 



47r[(6 + d ± 1)2 - 4d]r2 ' 



where h and d are non-negative constants. These solu- 
tions [2], as those of Pant and Sah [1], also suffer from 
the singularity problem since the matter density and the 
fluid pressure are infinite at the centre r = . This point 
for both solutions will be discussed later on. 

The total mass and charge contained within the 
sphere of radius a are, respectively. 



m 



{b + df + 2(6 - d) 



{h + d+iy - Ad 



a, 



(33) 



{h + df - 2d 



-,1/2 



a. 



(34) 



_{b + d+iy - Ad 

Now, vanishing of the charge in Eq. (34) implies that 

b + d = ±^/2d. (35) 

Again, by the use of the pure charge condition, i.e., the 
vacuum- fluid equation of state p+p = [22] in Eq. (32), 

one obtains 



b + d^O. 



(36) 



Thus Eqs. (35) and (36) imply both 6 = and d = 0. 
Therefore, substitutions of these values to Eqs. (30), (31) 
and (33) make all the physical quantities, including the 
gravitational mass, vanish and provide an electromag- 
netic mass model with the imperfect-fluid equation of 
state (i.e., the reduced Eq. (32) which now takes the 
form p + p = 0). 

However, when the constants b and d are identifled 

as 



b = [1 - c(n - l)^]/4c, 
d = [c{n + if - l]/4c. 



(37) 
(38) 



the solutions assumes the Pant and Sah form [1] which 
also provides an electromagnetic mass model, as we have 
verifled in Sec. 3.1. 

4. Conclusions 

In the present article, we have shown that some of the 
known solutions, related to static charged fluid spheres 
of Tolman- VI type, are of purely electromagnetic origin. 
In this regard, it is also to be mentioned that Ivanov's 
solutions [16] with a (A, n) classiflcation scheme are of 
this category. He has obtained a set of solutions by a 
Bessel function technique, which recover the Pant and 
Sah solution [1]. However, in this set of solutions, the 
density and charge function are singular at the centre, as 
stated earlier. It is also mentioned that the Tikekar solu- 
tions [2] suffer from the same singularity problem. One 
of the most undesirable features of general relativity is 
the occurrence of spacetime singularities where the laws 
of physics break down, and this is inevitable according 
to the Hawking-Penrose singularity theorems [33] . How- 
ever, following Herrera and Ponce de Leon [27], one can 
consider a sphere as being composed of (1) a central core 
of radius ro inside which all physical quantities are fi- 
nite, and (2) above the core, a self-similar fluid described 
by any of the solutions of Sections 3.1 and 3.2. Then 
these solutions can be matched with any suitable inte- 
rior Schwarzschild solutions across the surface r = tq 
where < ro < a. It is interesting to note that, in the 
context of cosmology, Trautman [34] suggested that in 
the Einstein-Cartan theory, where spin and torsion are 
taken into account inherently, a singularity is averted 
due to the action of torsion in a universe filled with 
spinning dust. 
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Regarding the stability of charged fluid spheres, Bon- 
nor [If] and also independently Do and Raychaudhuri 
[35] sliowed that a dust cloud of arbitrarily large mass 
and small radius can remain in equilibrium if it satis- 
fies the relation between the electric charge and matter 
energy densities by ct = ±p. However, Glazer [36], con- 
sidering radial pulsations, showed that the Bonnor [f f] 
model is electrically unstable. He [37] also explicitly es- 
tablished the effects of electric charge upon dynamical 
stability. According to Stettner [38], a fluid sphere of 
uniform density with a net surface charge is more stable 
than a neutral one. Whitman and Burch [39], for an 
arbitrary charge and mass distribution, showed that a 
charged analogue gives more stability. However, appli- 
cation of the pulsation equations to the charged Pant 
and Sah solution [f] gave an unsatisfactory result in 
connection with the boundary condition which is incom- 
patible with the densities and pressures [39] (a general 
technique for studying the stability is provided in the 
Appendix. 

We would like to point out that the solutions ex- 
pressed by Eqs. (20), (22) and (24) (except Eq. (26) for 
n = 2 of Pant and Sah [1]), do not obey the energy 
condition [33]. The reason behind this is related to 
the equation of state which, in the present case, takes 
the form p + p = Q due to application of the ansatz 
q{a)/a = 1/^/2. This equation provides two situations: 
either p > Q and hence p is negative or p > which 
means that p is negative. Thus, the strong energy con- 
ditions RijK^Kj > , where is a timelike vector, are 
violated in both cases. The first possibility mentioned 
above requires that the system should be under ten- 
sion and hence gravitationally repulsive in nature [21]. 
This case can be expressed in the form (?oo5ii = —1 and 
hence X = —jy , which is also equivalent to the charged de 
Sitter solution. The second option, viz., p <0 indicates 
a negative mass for the inside of the fluid sphere. This 
is possible in the case of Lorentz's extended electron of 
the size ~ 10~^^ cm where the spherically symmetric 
charged distribution of matter must contain some neg- 
ative mass density [40-43] . 

Appendix: Stability analysis 

To study the stability, one needs to match the interior 
solution to the exterior Reissner-Nordstrom black hole 
solution 



, , I 2m g \ , 9 2m q 



(A.1) 



at the junction interface S situated outside the event 
horizon, a > rh = m ± a/ttz^ — . To analyze the 
stability, it is required to use the extrinsic curvature, 
or second fundamental forms, associated with the two 
sides of the shell S as = —u^.^e'^.^e'^jy where 



are the unit normals to S and e'^'-j are the components 
of the holonomic basis vectors tangent to S. According to 
the Darmois-Israel formalism [44,45], one can write the 
Lanczos equations for the surface stress-energy tensors 
Sj at the junction interface S as [46-50] 



(A.2) 



where 5] = A\ag{—Q,pe,P(f,) is the surface energy ten- 
sor with Q the surface density and pe and the sur- 
face pressures; [Kij] = Kf- — K~- and K = [Kl\ = 
trace [Kij] . 

To analyze the dynamics of the junction shell, we 
permit the junction radius to become a function of 
proper time, a — > a(r). Now, taking into account 
Eq. (A.2), one can find 



1 

^ 47ra 
Pe = P4>=P 



2m 

a 

1 1 



mja -|- aa -|- 



87ra ^1 - 2m/a -h gVa2 -h d2 
{\-\-n){c-\-b?)-\-aa 



(A.3) 



(A.4) 



In what follows, the overdot means a derivatives with 

respect to r . 

The conservation identity 5*1. = — -I- 2^(p -|- q)\ 
yields the following relation: 



where 
Y = 



n 



Ana' 



Vc + d2. 



(A.5) 



(A.6) 



After a little bit of algebra, Eq. (A.3) gives the thin 
shell equation of motion 



d2 + v{a) = 0, 
where V{a) is the potential and defined as 

na) = l(/r+c)-4.2.V-i{L_£^, 



(A.7) 



(A.8) 



where /i = 1 — 2m/ a + q^/a^ ■ 

Linearizing around a static solution at a = oq , one 
can expand V{a) around ao to obtain 



V = V{ao)+V'{ao){a-ao) 

+ ^V"{ao)ia-aof+0[{a-ao)% 



(A.9) 



where the prime denotes a derivative with respect to a . 

Since we are considering linearization around a static 
solution at a = ao, we have V{ao) = and V'(ao) = 
0. Stable equilibrium configurations correspond to the 
condition y"(ao) > 0, i.e., V{a) has a local minimum 



Relativistic gravitational mass in Tolman- VI solution 



229 



at ao . Now, we define a parameter /3 , to be interpreted 
as the speed of sound [51], by the relation 



Following Eq. (A. 5), we get 



-^(jp + e} + y'-e" 



The solution is stable if 



Po < 



A-B+C-S-T+G-H 
N-L 



(A.IO) 



(A.11) 



(A.12) 



where /?„ = /3(a = oq) and A, B, C, S, T, G, H, N, 
L are given at a = oq as follows: 



H 



2-'^ 2 



6q 4m 
~ IF 



327r2£.3a2 

/2 



+ 



{h-cf 



1 (2m/a2 - 2gVa3)2 



T = 



G = 



(/i - 1 (/i - c)(-4m/a3 + e^Va^) 



327r2e2a2 
(/i-c)/{ 

(/i-c)(2m/a2-2gVa3) 



(^/A-^^)= 



g' g' 2 
— + — + - 

p a a 



1 

X - 



1 



/2 



1 f /2 



Ui-c? \'^q' 3g'' J_ 
167r2p2a2 [ ^ 2g^ ^ 2a'^ 



N 



3 2 1 / /2 
2a2 a2(^-^)Vv:^ 

+ A i (J^. 

4g^ (/i - cf 
a3 327r2£.3a2 



27ra5 



{h-cf 



2,2'K'^gg' 



= --Ay/h-V~c) 



(^/7T-^/3)3 



(A.13) 



whore (sec also (A. 8)) 

/i = l-2m/a + gV«^ 
/2 = l-2,m/a+2q'^/a'^. 

Thus if ao, w, c, m and g are specified quantities, 
then the stability of the configuration requires the above 
restriction on the parameter (3q . This means that there 
exists some part of the parameter space where the junc- 
tion location is stable. 
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